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Within the framework of generalized combinatorial approach, complexity is determined as a dis- 
order measure for hierarchical statistical ensembles related to Cayley trees possessing arbitrary 
branching and number of levels. With strengthening hierarchical coupling, the complexity is shown 
to increase monotonically to the limit value that grows with tree branching. In contrast to the 
temperature dependence of thermodynamic entropy, the complexity is reduced by the variance of 
hierarchical statistical ensemble if the branching exponent does not exceed the gold mean. Time 
dependencies are found for both the probability distribution over ensemble states and the related 
complexity. The latter is found explicitly for self-similar ensemble and generalized for arbitrary 
hierarchical trees. 

PACS numbers: 02.50.-r, 05.20.Gg, 05.40.-a 



I. INTRODUCTION 

Despite a daily appearance of hierarchy in society 
and comprehension of its role in physical, biological, 
economical and other complex systems [l]-[5] theory 
of hierarchically constrained statistical ensembles has 
been developed only at description of dynamics of spin 
glasses [BJ |7j. Formal basis of this theory is that 
states of hierarchically constrained objects are related 
to an ultrametric space whose geometrical image is 
the Cayley tree with nodes corresponding to statisti- 
cal (sub)ensembles |8J. The relaxation of hierarchical 
structures had been considered first [9j as a diffusion 
process on either uniformly or randomly multifurcating 
trees characterized by a diversity being a measure of 
the tree's complexity [TU]. Consequent study of the hi- 
erarchical ensembles has shown [11] their evolution was 
reduced to anomalous diffusion process in ultrametric 
space that arrives at a steady-state distribution over 
hierarchical levels, which is a Tsallis power law inher- 
ent in non-extensive systems |12j . Principle peculiarity 
of hierarchical systems consists in the splitting of each 
statistical ensemble residing on a given level into a num- 
ber of smaller subensembles with passage onto the lower 
level; subsequent descent down the tree results in more 
small subensembles (see Fig. [lj. 

From statistical point of view the set of above 
subensembles is characterized by the complexity, whose 
value determines disorder of the hierarchical coupling - 
in analogy with the entropy of thermodynamic systems. 
Formally, above quantities are defined equally, how- 
ever, their physical nature is absolutely different: if the 
entropy characterizes disorder in distribution of primi- 
tive structural units (for example, atoms), at the com- 
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statistical states 

FIG. 1: Characteristic form of the internal energy landscape 
(a) and related hierarchical tree (b) of complex system [13] . 



plexity definition their role is played by subensembles 
into which the whole statistical ensemble is subdivided. 
Consideration of self-similar hierarchically constrained 
ensembles within the framework of generalized combi- 
natorial approach shows [14] that with strengthening hi- 
erarchical coupling the complexity increases monotoni- 
cally to a certain boundary value; the latter decays with 
both the growth of the variance of this coupling and the 
reduction of a branching exponent of hierarchical tree. 
If the last of pointed dependencies is obvious (indeed, 
non-branching tree does not possess a complexity) , the 
reduction of the boundary complexity with hierarchi- 
cal ensemble scattering is meant to be abnormal, as 
in the conventional thermodynamic ensembles the en- 
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tropy always increases with temperature. Obviously, 
with branching change, promoting transformation of hi- 
erarchical ensemble into usual one, the pointed anomaly 
should be weakened and the complexity decrease with 
the ensemble scattering will pass into its increase. 

The present work is based on generalization of the 
statistical approach [T5] which allows, in particular, to 
establish the connection between the entropy of simple 
systems and the complexity of hierarchical ensembles. 
To avoid a misunderstanding, it is worth to note that 
the complexity introduced in Ref.[TU] characterizes a 
diversity of hierarchical tree's themselves, whereas we 
are aimed to consider the disorder of hierarchically con- 
strained ensembles in their statistical distribution over 
such trees. Along this line, in Section |TT] initial state- 
ments are given to define the distribution over states 
of self-similar statistical ensemble and the complexity 
of arbitrary hierarchical system. Main results are ob- 
tained in Section [rn| for the complexity definition within 
continuum approach, where the main contribution is 
shown to be given by the deepest levels of hierarchy. 
Accounting for the discrete character of distribution 
in Section |IV| shows that continuum approach derives 
semi-quantitative results already for several levels of hi- 
erarchy. Section [V] where the complexity definition is 
generalized for arbitrary hierarchical ensemble, is de- 
voted to discussion of the obtained results. Appendix 
contains details of generalized combinatorial approach 
being a basis of our consideration. 



II. MAIN STATEMENTS 

Generally, the behavior of a complex system is deter- 
mined by the cluster structure of the whole set of hier- 
archical levels, however, the property of self-similarity 
enables one to consider a typical cluster and the level 
number only. As a result, the system state is described 
by the probability pi = pi (t) to occupy the hierarchical 
level I whose distribution obeys the nonlinear Fokker- 
Planck equation [TTJ 



ToPi = 



d 
dl 



D 



dpi 
dl 



(1) 



Here, the level number I ^> 1 is supposed to be large 
to justify the use of the continuum approach, the over- 
point denotes the derivative with respect to the time t 
whose microscopic scale is r , and the diffusion coeffi- 
cient D represents the second moment of the intensity 
of transitions between microscopic states. Unlike the 
behavior of simple systems which is determined by the 
linear drift corresponding to the first moment of this 
intensity, the principle role in behavior of hierarchical 
ensembles plays the nonlinear term, fixed by both the 
factor e > and the exponent Q € [1,2]. At the ini- 
tial stage t < r d , r d = (e Q ~ 2 / D Q ~ l )l Q r , the contri- 
bution of diffusion is negligible and the characteristic 



scale of hierarchy increases with time according to the 
power law l c — Q 1 ^(t/To) 1 ^ , whereas the probabil- 
ity density Pl (t) = (Z/Qe) V(Q-i)(i/ To )-V(Q-i) decays 
hyperbolically with time. Transition to the diffusion 
regime r d ~ t <C t, t = 1 2 tq leads to transformation 
of the time dependence l c (t) into the usual root form 
l c = ^/2(i/Yo), while the probability density decays ac- 
cording to the same hyperbolic law. This law is ap- 
peared to be always inherent in hierarchical ensembles, 
being not only self-similar but also arbitrary ones [10] . 

With growth of the time to macroscopic values t^>r, 
the probability distribution takes the stationary Tsallis 
form \TB 



Pi = 



Po 



Po A 



2-Q 



(2) 



, A^D/e. 



The probability ^ increases monotonically with I de- 
crease, i.e., with growth of hierarchical cluster, reaching 
the maximum value po on the top level I = related to 
the whole system. Growth of the variance A = D/e ex- 
pands considerably the stationary distribution over hi- 
erarchical levels. Characteristically, in the limit A ^> 1 
the distribution Q differs slightly from exponential one 
on high levels I <C A 1 ^ 2- ^, however, with passing 
onto deeper levels the distribution tail becomes power- 
law. Study of possible types of hierarchical coupling 
[TJj] has shown the distribution ([2| is inherent in statis- 
tical systems related to a self-similar ultrametric space 
with fractal dimension d — (Q — l)/(2 — Q). 

The statistical theory of self-similar hierarchical en- 
sembles is based on generalization of both logarithmic 
and exponential functions given by expressions (A.ll, 
being type of deformed exponential distribution ([2J). As 
shown in Appendix, the additivity condition of the log- 
arithm function is provided by the related deformation 
rules for the multiplication/division operations, essen- 
tially changing the form of the usual multinomial coef- 
ficients which are the basis of combinatorial formalism 
of statistical physics [TS]. Making use of this formal- 
ism shows [14. that connection between complexities of 
the nearest hierarchical levels is expressed by equality 
(A.ll). In accordance with Fig{2]we suppose here that 
N n states of the top level are distributed over the m 
groups on the lower level (labeled by i), each of which 
contains nodes ij on the bottom level of a hierarchical 
cluster. If the group i possesses N% occupied states and 
there are Nij occupied states within each set ij, then the 
corresponding probabilities take the values Pi = Ni/N, 
Pij = Nij/N, where N is the total number of statisti- 
cal states of hierarchical ensemble. Note that N is not 
reduced to value N <§; N related to the given cluster. 
Because of the obvious condition Ni = Y^Tlx^iji tnc 
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FIG. 2: a) Numbers of states on nearest levels of typical 
hierarchical cluster; b) trivial hierarchical tree; c) parame- 
terization of nodes on nearest levels. 



above probabilities are connected by the equality 

mi 

Pi = E^- 



(3) 



As a result, expression (A.ll), accompanied with Eqs 



(A.8I (A.9), takes the form 



Cq (pu, . . . ,p mmm ) = Cq (pi, . . . ,p rn ) + 



E*?<* 



pa 
Pi ' 



pi 



(4) 



Here, the last factor represents relative complexity de- 
fined by the Tsallis expression 



C c 



Pa 



Pi 



(Q-i)- x E 



PHj_ 

p% 



Pi J 



Q 



(5) 



Utilizing this formula we express the connection be- 
tween the complexities of nearest hierarchical levels in 
terms of the related probability distributions over sta- 
tistical states: 



Cqipu,... i Pmm m 



Q 



1 m mi 

^tEE^ 1 



i=i j=i 



(6) 



Above consideration concerns a typical hierarchical 
cluster depicted in Fig(2} To consider the whole sta- 
tistical ensemble consisting of such clusters, let us sup- 
pose that statistical states are distributed over micro- 
canonical (sub)ensembles with probabilities and corre- 
sponding complexities defined by the level number I 
only: {py} p l+1 , {p. t } =$> pi; Cq (p n , . . . ,p mmm ) => 
C(l + 1), Cq (pi, . . . ,p m ) C(l). As a result, relation 
(l6| takes the simplest form 



C(l + l)-C(l) 



M ( Q-l Q-l\ in\ 

-Pi+i [Pi -Pi +1 ) ■ (7) 



Q-l 



For self-similar ensembles node numbers are distributed 
over levels I = 0, 1, . . . , n according to the power law [TS] 



Mi = (l+ l) a 



(8) 



with an exponent a > 1. In the limit / S> 1, the condi- 
tion ([3| takes the form 



Mipi = M l+1 pi +1 . 



(9) 



Then, Eqs. (|8]), ([9| accompanied with asymptotic be- 
havior of distribution ([2| show the self-similarity con- 
dition is fulfilled for hierarchical trees related to the 
branching exponent 



(10) 



Q-i 



Stationary probability distributions over states of the 
nearest levels are known to be connected by the recur- 
rence relation [11] 



Pi+i - Pi 



-pf/A, 1 = 0,1, 



(11) 



In the limit I ^> 1, this relation is reduced to the steady- 
state Fokker-Planck equation 0. 



III. CONTINUUM APPROACH 

Assuming that total number of levels is n — > oo, let us 
find the complexity of hierarchical ensemble where the 
main contribution is given by deep levels with / ^> 1. 
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Here, the statistical state probabilities of nearest levels 
are connected by relations 



Q-l 
Pi+i 



Pi 



Q-i 



d Q-x 

~ dA 
(Q - i)p^"i 



,dp 



i+i 



dl 



Q-l 2(Q-1) 



(12) 



where the last equality is written with accounting for 
equation ( 11 I taken in the continuum limit. Then, sub- 
stitution of Eq.(Jl2) into Eq.Q gives 



C(l + 1)-C(l) 



M, 
A 



l 2Q-1 
-Pl+1 ■ 



(13) 



Rewriting C(l + 1) — C(l) in the continuum limit as 
dC/dl and taking into account Eqs. Q allow one 
to reduce the difference equation ( [13] ) to the following 
differential equation, 



dC 
~dl 



I" 
A 



-(Q-l) 



2Q-1 



(14) 



where I + 1 is replaced with I . 

In the limiting case a — > 1, solution of Eq.(14l leads 
to the dependence 



C(n) = ^ 



Pn\ 

PoJ 



(15) 



Here, the complexity appears to be the monotonically 
increasing function of the level number n. It has the 
initial value C(0) = and the maximum value 



C(oo) 



(2-Q)^ 



A -3=1 

-A 2 -Q. 



(16) 



decaying with growth of both the non-extensivity pa- 
rameter Q and the variance A. Remarkably, the com- 
plexity obtained here behaves similarly to that found in 
our recent work |14j . However, the analogous expres- 
sions (16)-(18) appearing in Ref. |14| are much more 



complicated than the equalities ( 14 1-( 16 ) derived in this 
section. It is accounted for that continuum approach 
has been already used in Ref. [14] at estimation of the 
relative entropy whereas in the present study we 
pass to continuum limit at a later stage, while defining 
probability difference (12) to estimate the complexity 
difference Q. At the same time, no approximations 
have been imposed to express the entropy The 
pointed distinction resulted in the exponent of the dis- 
tribution function derived in Ref. [T3] (see the r.h.s. of 



Eq.(16) therein) to exceed the exponent value in Eq.( 14) 
of this section by the magnitude Q — 1 < 1. This fact, 
however, does not play a principle role within our ap- 
proach. 



At arbitrary values a of the branching exponent, so- 
lution of equation ( 14 1 is expressed by means of the 



hypergeometric function: 
CW^^(^,l + fl ;2 + fl ;_ 

_ Q-i q-i 



A 



-Po 



(17) 



The analysis of r.h.s. of Eq.(14l shows that its solution 
holds finite values in the limit n — > oo if the exponent 
a does not exceed the maximum value 



Q 



Q- l 



(18) 



As shown above, the self-similarity condition leads to 
the lower value ( 10 1 for a, at which odd arguments of 



the hypergeometric function coincide and one arrives at 
the binomial form [17] 



F 



2Q - 1 Q 2Q - 1 



-1 'Q-l' 



- 1 



(1 + vri) Q - 1 



(19) 

As a result, the complexity (17) is expressed by the 
simple equation [20 



C(n) = 



Q 



pz nQ - 



(20) 



As for the degenerate tree (a — » 1), the complexity in- 
creases monotonically with the number of hierarchical 
levels to the maximum value 



C(oo) 



(2 



, Q-l 



Q(Q-i)^ 



Q+IQ-I , n 
A (Q-i)(2-Q) \&+ 



Q) 



(21) 



where roots Q± = (1 ± v5)/2 represent the gold mean. 
With growth of the non-extensivity parameter Q, the 
maximum complexity decreases monotonically from in- 
finite value at Q = 1 to zero at Q = 2. However, un- 
like Eq.(16l, the complexity dependence on the vari- 
ance A becomes non-monotonic: at values of the non- 
extensivity parameter limited to the top magnitude 
Q + ~ 1.618, the maximum complexity increases with A 
growth, whereas at Q > Q+ it decreases. With account- 



ing for the self- similarity condition (10), this means 



that behavior inherent in simple statistical systems is 
realized at the branching exponent of hierarchical tree 
exceeding gold mean a + = (Vb + l)/2 ~ 1.618; on 
the other hand, the complexity decay with the variance 
growth, what is characteristic for complex systems, is 
shown to appear at values of the branching exponent 
within the domain 1 < a < 1.618. 
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IV. DISCRETE HIERARCHICAL LEVELS 

For hierarchical systems possessing a finite number 
of levels n > 1, the above results of the continuum ap- 
proach lose accuracy, and we need to take into account 
the discretization of hierarchical statistical ensemble. 
Probability distribution of such ensemble is defined by 
the set of equations (111, whose number n is less by one 
than number of probabilities pi, I = 0, 1, . . . , n. The en- 
semble description is reached by adding to the system 



(11 1 the normalization condition 



(22) 



Numerical solution of Eqs. ( [T7| ), p2) arrives at the 
probability distributions over hierarchical levels with 
different variances, as shown in Fig [3] (for convenience 
we use both half- logarithmic axes, where exponential 
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FIG. 3: Probability distributions over hierarchical levels 
within representations of both half-logarithmic (a) and log- 
log (b) axes at Q = 1.5 (curves 1-6 relate to values A = 
0.9, 1.25, 2.0, 3.0, 4.0, 6.0, respectively). 

dependencies are straightened, and logarithmic axes, 
representing power laws with straight lines). As Fig|3]i 
shows, distributions of statistical states over hierar- 
chical levels are far from exponential dependencies at 
all values of the variance A. The pronounced curve 



straightening in FigjSJ) indicates the transition into 
power-law regime in the limit I ~ * oo. The latter is 
confirmed by comparing the distribution function (12) 
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FIG. 4: Comparison of the distribution function |2) found 
within continuum approach (solid lines) with solution of the 
difference equations and the normalization condition 
(22) (points) at Q = 1.5, A = 1 (a) and Q = 1.5, A = 5 (b). 



derived in the continuum limit with the solution of the 



system of difference equation (111 subject to the nor- 
malization condition (22 1. As illustrated in Figjl] both 



approaches produce virtually the same results for the 
deep levels I ^> 1, however, the discrepancy at lower 
level numbers becomes bigger as the variance A in- 
creases. 

The most precise representation of accuracy of the 
continuum approach is reached if one expresses the dis- 
tribution over hierarchical levels in the modified Tsallis 
form 



Pi 



Po 



Po + A ' 



2-q 



(23) 



< I < n 



with the exponent ai taking the value ai ~ Q — 1 in the 
limit I — > oo. As Fig(5] shows, the continuum approach 
is improved by both the tendency of the non-extensivity 
parameter Q to the limit value Q — 2 and the unlimited 
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FIG. 5: Dependencies of the exponent of distribution ([23f 
on the level number at A = 1 (a), A = 2 (b) and Q = 1.5 
(c) (numbers near curves point to values Q in panels (a), 
(b) and A in (c)). 



growth of the variance A. In particular, Q — > 2 means 
that branching exponent ( 10 1 tends to the limit a = 1 



related to the degenerate hierarchical tree, where the 
statistical ensemble is defined by hyperbolic distribu- 
tion with the exponent ai ~ Q — 1 — > 1. Respectively, 
growth of the variance increases the disorder of statis- 
tical ensemble dramatically and the discrete nature of 
the ensemble becomes inessential. 

As in stationary case, evolution of the probability dis- 
tribution over hierarchical levels is determined by the 



condition ( 22 ) . The number of the deepest level is de- 
fined by the diffusion law n — W2(t/TQ) (for the sake of 
simplicity, we accept the condition Ta/r <C 1 for charac- 
teristic time scales, when the contribution of abnormal 
drift is negligible). Corresponding time dependencies 
of the probability distribution functions are depicted in 
Fig|6] It is demonstrated how diffusion process over 




FIG. 6: Evolution of the probability distribution over hier- 
archical levels at Q = 1.5 and A = 3 (curves 1-4 relate 
to the time moments t/ro = 50,1.25 ■ 10 3 ,5 ■ 10 3 ,5 • 10 5 , 
respectively). 

levels of the hierarchical tree provides relaxation to the 
stationary distribution shown in Fig|4] 

Let us proceed now with definition of the complexity. 
Taking successfully values I = 0, 1, . . . , n for the number 
of levels in Eq.([7]), with accounting for Eqs. (|8]), (10 1, 
one obtains 



C(n) = 



1 



1 ^ 

i=i 



*/i/(Q-i) 



Pi [P^ 



■Pi 



Q-l 



(24) 



As shown in Fig|7] with strengthening hierarchical cou- 
pling the complexity grows monotonically to the maxi- 
mum value which decays to zero as the non-extensivity 
parameter Q tends to the upper border Q — 2. What 
about the complexity dependence on the variance A, it 
is much more complicated: at small values of Q the 
A-growth promotes the complexity-increase, whereas 
at large Q the complexity decreases with the vari- 
ance. Fig(8] demonstrates this peculiarity representing 
the complexity as a function of variance A for a large 
number of levels n = 500 and different values of Q. It is 
seen that complexity increases monotonically with A at 
small Q, decays at large Q and varies non-monotonically 
at intermediate Q. 

According to consideration of the previous Section, 
such complexity behavior is also captured by the rela- 
tion (21 1. The latter, in particular, implies that the 



system of equations (11 1 along with the normalization 



maximum complexity as a function of variance exhibits 
the behavior inherent in simple statistical systems at 
values of non-extensivity parameter limited from above 
by the gold mean Q + ~ 1.618. On the other hand, 
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FIG. 7: Dependence of the complexity on the number of 
hierarchical levels at Q — 1.3 (a), Q = 1.5 (b) and Q — 1.7 
(c) (numbers near curves point to the A values). 



at Q > Q+ the variance-growth reduces the complex- 
ity of hierarchical ensemble. Qualitatively, such com- 
plexity dependence can be perceived while consider- 
ing a trivial hierarchical tree (see Fig[2]3) in line with 
the following reasonings. As the top node splits into 
to > 1 parallel branches, each of which comprises of 
77 > 1 sub-nodes, the total number of nodes equals 
N = 1 + ran; at that, the top node relates to the 
probability po = m/N, whereas each of the remaining 
sub-nodes is associated with the probability p — 1/N. 
Then, the complexity definition (24 1 gives the value 

Cx = (Q - l) _1 i>(Po _1 - for each of the tree 
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FIG. 8: Dependence of complexity on variance for levels 
n = 500 at Q = 1.50 (a), Q = 1.65 (b) and Q = 1.70 (c). 



branches and the total complexity C = mC\ takes the 
value 



C=(Q-1) 



-il 



-(Q-l) 



(77 + 777 1 ) c 



(25) 



As one would expect, the latter grows with increase of 
the branching 777 and reduction of the length of branches 
77. Since for a stochastic tree the 77- growth is attributed 
to the increase in disorder of states over levels, the es- 
timation ( |25[ ) provides an explanation for an abnormal 
decay of the complexity with growth of the variance of 
hierarchical ensembles - it is caused by the increase in 
number of subensembles which are connected with each 
other without hierarchical constraining. 
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Characteristic peculiarity of the curves depicted in 
Fig. [8] is their break at small variance, where the 
complexity (24) becomes ill-defined. Such behavior is 
caused by the functional form of the distribution 
implying that the decrease in variance A leads to un- 
limited growth of the probability to occupy the top level 
Po oc A _1 /( 2_< 2). Therefore, at small A the maximum 
probability p — 1 is reached, so the whole hierarchical 
system accumulates on the top level and the definition 
of complexity loses sense. 



Comparison between the complexities (24 1 and (20) 



defined within discrete and continuum approaches, re- 
spectively, is shown in Fig(9j It is seen that both ap- 




FIG. 9: Comparison of complexities (24 1 determined within 



both discrete (solid lines) and continuum (dashed lines) ap- 
proaches: a) Q = 1.5 (curves 1-3 relate to A = 1,2,3, 
respectively); b) Q — 1.7 (curves 1-3 relate to A = 
0.7,0.9,1.2, respectively). 

proaches result in the identical complexity dependence 
on the number of hierarchical levels n. However, with 
growth of the non-extensivity parameter Q the contin- 
uum approach loses the accuracy. Such behavior is ra- 
tionalized by the reduction of the tree branching with 
Q-growth, promoting enhanced contribution of top lev- 
els where the continuum approach is inapplicable. 

Time dependence of the complexity C(t) is defined 
by evolution of the probability distribution pi(t) over 
hierarchical levels. Fig(6] shows that in the course of 
time the level variance is strengthened, so, accordingly, 



the dependence C(t) should have the same form as the 
complexity variation with the A variance growth for 
stationary case (see Fig|8]). Indeed, as indicated in 
Fig |10| complexity increases with time the faster, the 
greater the distribution scattering, provided the non- 
extensivity parameter does not exceed the gold mean 
Q + ~ 1.618 [13]. However, at Q > Q + the system evo- 
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FIG. 10: Time dependencies of the complexity at Q = 1.5 
(a) and Q = 1.7 (b); curves 1-3 relate to values A = 1, 2, 3 
(a) and A = 0.7, 0.9, 1.2 (b), respectively. 

lution becomes abnormal: a short-term increase of the 
complexity from finite values is followed by its long- 
continued decay to the values decreasing with the A 
variance growth. 



DISCUSSION 



Above consideration is based on the definition of mi- 
crocanonical ensemble, where the probabilities Pi 1 ...i l to 
occupy statistical states i\ . . . ii on Zth hierarchical level 
take equal values pi ~ po/Mi defined by the condition 
([£]). According to Eq.Q, these probabilities decay with 
Z-growth due to increase in the number of states Mi . As 
a result, the recurrence relation ([6| is reduced to sim- 
ple equality (|7|, whose iteration gives the expression 
(24 1 for the complexity of the self-similar hierarchical 



ensemble. 
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In general case of arbitrary distribution over nodes 
of the hierarchical tree, not possessing the property of 
self-similarity, the difference equation (G) leads to the 
complicated expression 



1 " 

^ 1=1 

m m il raii* 2 ....,_i 

^1 = ^2^2 " E Pii... it-iii x 
— 1 22 — 1 it=l 



fined as 



(26) 



( P 9- x 



"in... 



where the probability Pi 1 ...i l _ 1 is reduced to the value 
po at £ = 1. In contrast to a simple two-level tree de- 
picted in Fig(2j it is meant here that on the level I of 
the hierarchical tree comprised of n > 1 total number 

of levels the set of states i\ . . . i\ . . . ?'/-i2 

i\ . . . ii-imi 1 i 2 ...ii_ 1 forms a cluster related to the node 
i\ . . . ii-i on the upper level l—l (see Fig{2j;). Therefore, 
to calculate the complexity ( 26 1 , one needs first to sum 
over the nodes i\ . . . i n -\i n of the bottom level n, which 
belong to the cluster related to the node i\ . . . i n -i of 
the upper level n — 1. Then, the summation over nodes 
of the level n belonging to all remaining clusters is car- 
ried out and the same procedure is repeated for each 
subsequent hierarchical level I < n [21] . 



Expression ( 26 ) is a basis for the numerical defini- 



tion of the complexity of arbitrary hierarchical ensem- 
ble (for example, in case of the complex defect structure 
of solids subject to intensive external influence type of 
strong plastic deformation or rigid radiation treatment). 
Unlike the amorphous systems, the number of structure 
levels of a real crystal is rather not large. Particularly, 
among different spatial scales in the latter type of sys- 
tems, it is accepted to distinguish between several basic 
levels of consideration [13]. The finest level, the mi- 
croscopic one, is related to homogeneous distribution of 
point defects, dislocation and disclination. Next, the 
coarser mesoscopic level corresponds to homogeneously 
distributed structure cells and fragments. Finally, the 
macroscopic objects such as homogeneously distributed 
grains and texture components are considered on a 
macroscopic level. To define the complexity of a real 
structure, one needs first to distribute the whole en- 
semble of defects over hierarchical levels I = 0, 1, . . . , n, 
then to calculate on each of them a number of defects 
belonging to the cluster i\ . . . i% . ■ ■ «/-i2, 



. . . , %\ . . . %\— lTTii 1 i 2 ,..ii_i , 

probability 



Pii.. Al- 



and, finally, to attribute the 



N 



(27) 



to the node i\ . . . of the upper level l — l. At that, 
the total number of defects residing on all levels is de- 



ii—li-2 — l in — 1 



where 



_j is the number of possible states i\ 



within the cluster related to the node (at 
n = 1 one has mj 1 ...i n _ 1 = m) |22j . Substitution of 
the obtained set of probabilities {pi 1 ...ii }, I — 0, 1, . . . , n 
into Eq.( |26[ ) gives the complexity of hierarchically con- 
strained defect structure of a solid. Obviously, this com- 
plexity determines such phcnomcnological quantities as 
strength and plasticity of solid. Remarkably, such def- 
inition of the structural complexity cannot be reached 
solely by means of experimental methods (for example, 
electron microscopy), but requires a subsequent com- 
puter processing according to the above algorithm. 

Realization of the described program, however, rep- 
resents rather a challenge. Therefore, in the preced- 
ing sections we limited ourselves to considering a self- 
similar hierarchical tree characterized by the level num- 
ber n and a typical cluster with the branching exponent 
a. Studying the various hierarchical trees has revealed 
three their types [15]: (1) the degenerate tree possessing 
only one branching node per level, so the total number 
of nodes increases linearly with n; (2) the regular tree, 
where on each level all nodes branch equally, and the 
total node number increases exponentially with n; and 
(3) the self-similar irregular tree characterized by the 
power law ^ with exponent a > 1 for the Z-dependence 
of the number of nodes. In the first two cases the 
probabilities of state occupation on various hierarchi- 
cal levels vary logarithmically slow and exponentially 
fast, respectively, while in the last case it obeys the 
power-law dependence Q inherent in the self-similar 
systems. According to Eq. ([To]) , the non-extensivity pa- 
rameter Q — (a + l)/a is determined by the branching 
exponent a: regular (a = oo) and degenerate (a = 1) 
trees are characterized by limiting values Q — 1 and 
Q — 2, respectively, whereas the power-law distribution 
([2]) with exponent 1 < Q < 2 takes place at 00 > a > 1 . 

As pointed in Introduction, the complexity of differ- 
ent hierarchical trees has been first considered in works 
[9], [HJ|. Within our notation, these trees have been 
characterized with the silhouette 



si 



In 



Mi 



Ml 



(29) 



that determines the logarithmic growth rate of the node 
number Mi with increase of the hierarchical level num- 
ber I. A peculiarity of the approach used in Rcf.[9 
is that it is based on the consideration of the regular 
tree, whose nodes are multifurcated with the constant 
branching index b > 1. In this cas e, o ne has the de- 
pendence Mi = b l and the definition ( 29 1 gives a simple 
relation s = In b. However, for the main object of our 
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consideration, a self-similar tree, Eq. ([8| for the num- 
ber node distribution has to be used, and the resulting 
silhouette s/ = In (1 + l/l) a ~ a/ 1 depends significantly 
on the level number even in the continuum limit I 1 . 
It is clear that such /-dependence is caused by the def- 



inition ( 29 1 introduced for trees related to the regular 



ones, whereas self-similar irregular trees represent their 
antipode. Obviously, the silhouette of a self-similar tree 
may be defined with the equalities 



M xl - Mi 
(A - 1)M, 



1 



(30) 



This definition is based on use of the Jackson derivative 



which determines the variation rate of the function ( 10 ) 



with respect to the dilatation A > 1 and reduces to the 
usual derivative in the limit A — > 1 [18 . According to 



Eq.(29l, the silhouette of self-similar tree is the A-basic 



number [a]\, whose value equals the exponent a at the 
limit dilatation value A = 1 and grows as A a_1 at A ^ 1. 
Obviously, the present study of self-similar hierarchical 
ensembles relates to the first case, A = 1. 

In summary, our study indicates that formation of hi- 
erarchical coupling promotes a fast complexity growth 
to a certain maximum value. With strengthening of 
the branching parameter a of hierarchical tree the max- 
imum complexity grows monotonically from zero at 
a = 1 to infinity at a — > oo. Investigating the variance- 
dependence of the complexity has revealed a behavior 
inherent in simple systems at the values of branching 
exponent exceeding the gold mean a+ = 1.618. On the 
other hand, the complexity-decrease with the variance- 
growth, being a characteristic of complex systems, has 
been observed for the range 1 < a < 1.618 of the 
branching exponent. 



Appendix 



As known, the non-extensive statistical mechanics is based on the following definitions for logarithmic and 
exponential functions [12 . 

ln q (x) := X 9 ~ 1 , exp (a?) := [1 + (1 - q)x]^ ; [y]+ := max(0, y), q < 1, (A.l) 
1 — q H 

which are reduced to the usual functions in the limit q — > 1 . Introducing the g-deformed multiplication and division 
operations for the positive values x, y as follows, 



x ®, y = [x l ~ q + y 1 -" -l]]- q , i0,!/= [x l ~ q - y^ q + l] ^ J 



x,y > 0, 



it is easy to verify that they satisfy the usual properties 

ln g (a; <& q x) — ln q x + ln q y, ln q (x <Z) q x) = \n q x - ln 9 y; 

exp q (x) ® g exp ? (y) = exp q (x + y), exp q (x) (Z) q exp q (y) = exp q (x - y) 



(A.2) 



(A.3) 



of the logarithmic and exponential functions. 

Within combinatorial approach |T5], the g-deformed statistics is reduced to consideration of the generalized 



factorial N\ q := 1 ® q ■ ■ ■ ® q N and the corresponding logarithm, 



\n q {N\ q ) 



i-q 



(A.4) 



In the thermodynamic limit N — > oo, a sum in the above equation is replaced by an integral, and one gets 



Defining g-deformed multinomial coefficient as 

N 

JVi . . . N k/ 



(M q )0 q [(JVi!,)®,"-®, (JVfe!,)] 



with a set of integers Ni subject to the condition N — X)"=i we f^ n< ^ 



N 

Nx ... N k 



N 



N k 



E^ 9 -E^---E^+ 1 



(A.6) 



(A.7) 



J + 
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Similarly to Eq.(A.5l, the logarithm 



N 



ln„ 



N 



2-9 



k 



is calculated to construct a Tsallis entropy 

Cq(pi, ...,Pn 



-d(N)+ ECi(Ni), 



N V W r^ 1 - 1 



9 = 2 



(A. 



(A.9) 



related to the physical parameter Q = 2 — q > 1. 

The above formalism can be easily generalized to study hierarchical systems [19]. To this end, let us consider a 
structure of the hierarchical ensemble comprised of N states. These states are distributed over subensembles i — 
1,.. . , m, each containing Ni states. In turn, every of these subensemble further splits in m, smaller subensembles 
ij, each possessing states. As the ensemble states are distributed, the relations X^i Nij = Ni, TlT=i Ni = N 
hold. Then, expression (A. 6) for the multinomial coefficients takes a generalized form, 

N 



N 

Nu... N mm 



N...N„ 



Ni 

Nn ■ ■ ■ N lm] 



whose g-logarithm is 



ln„ 



N 

Nn... N r , 



ln„ 



N 

Ni...N v 



9 i=l 



N a ... N mt 



(A.10) 



(A.ll) 



Applying the estimation (A. 8) to the last formula arrives at the connection Q between the complexities of the 
nearest hierarchical levels. 
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